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This  paper  develops  a  mathematical  theory  for  solving  determin¬ 
istic,  Lanchester-type,  ‘square-law’  attrition  equations  for  combat  be¬ 
tween  two  homogeneous  forces  with  temporal  variations  in  fire  effec¬ 
tivenesses  (as  expressed  by  the  Lanchester  attrition-rate  coefficients). 
It  gives  a  general  form  for  expressing  the  solution  of  such  variable- 
coefficient  combat  attrition  equations  in  terms  of  Lanchester  func¬ 
tions,  which  are  introduced  here  and  can  be  readily  tabulated. 
Different  Lanchester  functions  arise  from  different  mathematical 
forms  for  the  attrition-rate  coefficients.  We  give  results  for  two  such 
forms:  (1)  effectiveness  of  each  side's  fire  proportional  to  a  power 
of  time,  and  (2)  effectiveness  of  each  side's  fire  linear  with  time  but 
with  a  nonconstant  ratio  of  attrition-rate  coefficients.  Previous  re¬ 
sults  in  the  literature  for  a  nonconstant  ratio  of  these  attrition-rate 
coefficients  only  took  a  convenient  form  under  rather  restrictive  con¬ 
ditions. 


|~^ETERMINISTIC  Lanchester-type  equations  of  warfare  (see  references 
38,  44)  are  of  value  for  identifying  trends  in  weapon  system  analysis  or 
force  structuring  studies  because  of  their  computational  convenience,  even 
though  combat  between  two  opposing  military  forces  is  indeed  a  complex 
random  process  (see  Note  1).  In  this  paper  we  present  a  mathematical 
theory  (including  a  new  standard  form)  for  the  solution  of  variable-coeffi¬ 
cient  Lanchester-type  equations  of  modern  warfare  for  combat  between  two 
homogeneous  forces.  After  some  preliminaries,  we  give  a  more  precise 
statement  of  this  paper’s  purpose  at  the  end  of  Section  3. 

First,  we  review  Lanchester’s  classic  mathematical  model  of  combat  be¬ 
tween  two  homogeneous  forces  and  its  extension  to  cases  of  time-varying 
fire  effectiveness.  Then  we  discuss  previous  work  on  developing  analytic 
solutions  to  variable-coefficient  Lanchester-type  equations  and  explain  how 
we  extend  these  results.  Next,  we  develop  a  mathematical  theory  for  solv¬ 
ing  Lanchester-type  equations  for  a  ‘square-law’  attrition  process  for  com¬ 
bat  between  two  homogeneous  forces.  These  general  results  are  then  ap- 
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plied  to  cases  of  some  particular  attrition-rate  coefficients  of  interest,  and 
new  mathematical  functions,  Lanchester  functions,  that  could  be  tabulated 
in  the  future  are  proposed.  Some  numerical  examples  are  given.  Finally, 
we  discuss  the  development  of  tabulations  of  values  of  Lanchester  func¬ 
tions. 


1.  LANCHESTER’S  CLASSIC  FORMULATION 

F.  W.  Lanchester  was  an  English  aeronautical  engineer  who  lived  from 
1868  to  1946  (see  reference  25).  In  1914  he  hypothesized  in  reference  23 
that  combat  between  two  opposing  forces  could  be  modelled  by  (see  Note  2) 

dx/dt=—ay  with  x(t  =  O)  =  x0, 

dy/dt=  —bx  with  y(t  =  0)=y0, 

where  t  =  0  denotes  the  time  at  which  the  battle  begins,  x(t)  and  y(t)  de¬ 
note  the  numbers  of  X  and  Y  at  time  t,  and  a  and  b  are  nonnegative  con¬ 
stants,  which  are  today  called  Lanchester  attrition-rate  coefficients.  These 
coefficients  represent  the  effectiveness  of  each  side’s  fire  (i.e. ,  its  firepower). 
Lanchester  considered  his  now  classic  combat  formulation  ( 1 )  in  order  to 
provide  insight  into  the  dynamics  of  combat  under  “modern  conditions” 
and  justify  the  principle  of  concentration  (see  Note  3).  The  above  equa¬ 
tions  are  valid  only  for  x,  y>  0.  The  first,  for  example,  becomes  dx/dt  =  0 
for  x-0. 

We  shall  refer  to  the  equations  (1)  as  Lanchester’ s  equations  of  modern 
warfare.  Two  sets  of  circumstances  under  which  these  equations  have  been 
hypothesized  to  apply  are : 

(a)  both  sides  use  aimed  fire  and  target  acquisition  times  are  constant, 
independent  of  the  force  levels  (a  special  case  is  when  they  are 
negligible), 

(b)  both  sides  use  area  fire  and  a  constant  density  defense. 

A  more  complete  discussion  of  these  hypotheses  is  to  be  found  in  references 
11  and  44.  Other  factors  may  be  included  in  the  equations  and  other  dif¬ 
ferential  equation  models  of  combat  attrition  may  be  referred  to  as  Lan- 
chester-type  equations,  but  we  will  not  consider  these  here  (see  references 
15  and  38). 

From  ( 1 )  Lanchester  deduced  his  classic  square  law 

b[x02-x\t)]  =  a[y02-y2(t)],  (2) 

which  has  the  important  implication  that  a  side  can  significantly  reduce  its 
casualties  by  initially  committing  more  forces  to  battle.  We  also  have  the 
well  known  result  that  the  time  history  of,  for  example,  the  X  force  level  is 
given  by 

x(t)  =  x0cosh(ab)~st  —  y0(a/b)ismh(ab)H.  (3) 
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One  result  of  this  paper  is  a  generalization  of  (3)  to  the  general  case  of  vari¬ 
able  attrition-rate  coefficients. 

2.  VARIABLE  ATTRITION-RATE  COEFFICIENTS 

The  pioneering  work  of  Bonder 16,7  401  on  methodology  for  the  evaluation 
of  military  systems  (particularly  mobile  systems,  such  as  tanks,  mechanized 
infantry  combat  vehicles,  etc.)  provides  a  motivation  for  interest  in 
variable-coefficient,  deterministic,  Lanchester-type  combat  formulations. 
Bonder161  has  pointed  out  that  (at  least  for  the  case  of  mobile  weapon  sys¬ 
tems)  the  validity  of  the  assumption  of  constant  attrition -rate  coefficients 
is  open  to  question.  Let  us,  therefore,  consider 

dx/dt—  —a{t)y  with  x(t  =  0)  =  x0, 

(4) 

dy/dt=  —b(t)x  with  y(t  =  Q)  =  y<>, 

where  a(t)  and  b(t)  denote  time-dependent  Lanchester  attrition-rate  co¬ 
efficients.  Recently,  Taylor  and  Parry  have  shown  by  consideration  of 
the  force-ratio  equation  that  a  generalization  of  (2)  holds  even  for  the  gen¬ 
eral  case  of  variable  attrition-rate  coefficients  [see  equation  (4)  of  reference 
41].  Accordingly,  we  will  also  refer  to  (4)  as  the  equations  for  a  ‘square- 
law’  attrition  process  (see  also  references  37  and  40). 

Two  significant  developments  in  the  Lanchester  theory  of  combat  during 
the  1960’s  were  (a)  the  development  of  methodology  for  the  prediction  of 
Lanchester  attrition-rate  coefficients  from  weapon-system  performance 
data  by  Bonder16 '81  and  others  in  references  3  and  10,  and  (b)  Clark’s1121 
development  of  methodology  for  the  (maximum  likelihood)  estimation  of 
such  coefficients  from  Monte  Carlo  simulation  output  data.  Both  of 
these  developments  and  others130,311  have  facilitated  the  application  for  de¬ 
fense  planning  studies  of  models  such  as  (4)  and  its  generalization  to  com¬ 
bat  between  heterogeneous  forces  (see  reference  10). 

The  Lanchester  attrition-rate  coefficients  in  (4)  depend  on  a  number  of 
variables,  such  as  firing  doctrine,  firing  rate,  rate  of  target  acquisition,  force 
separation,  tactical  posture  of  targets  (see  reference  5  or  pp.,  18-26  and 
pp.  81-114  of  reference  10  for  attrition-rate  coefficient  prediction  method¬ 
ology).  Bonder16,101  has  considered  a  number  of  forms  for  the  range  de¬ 
pendency  of  attrition-rate  coefficients  based  on  examination  of  data  for 
some  representative  weapon  systems.  Motivated  by  this  work,  we  will 
consider  the  following  coefficients : 

(I)  a(t)=ka(t+C)m  and  b(t)  =kb(t+C)n  with  C^O,  (5) 

(II)  a{t)  —ka(t+C)  and  b(t)  =kb(t+C+A)  with  A,CA 0,  (6) 

where  the  constant  A  allows  the  ratio  of  attrition-rate  coefficients  that  both 
vary  linearly  with  time  to  be  nonconstant  and  in  (6)  the  constant  C  allows 
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the  fire  effectiveness  of  each  of  the  two  systems  to  be  nonzero  at  the  begin¬ 
ning  of  battle  (i.e.,  <  =  0).  In  (5)  the  modelling  role  of  C  depends  on  the 
values  of  m  and  n;  for  example,  if  m  and  n>0,  then  C  allows  the  initial  fire 
effectiveness  of  each  system  to  be  nonzero.  Some  situations  that  may  be 
modelled  with  these  coefficients  are  discussed  below.  [Additionally,  for 
m  and  n  negative  (for  m,n^  —  1  we  require  C>0),  equations  (4)  and  (5) 
may  be  considered  to  model,  for  example,  an  infantry  ‘fire  fight’  in  which 
the  combatants  take  cover  so  that  the  effectiveness  of  fires  decreases  with 
time.]  For  a  battle  in  which,  for  example,  force  separation  or  range  be¬ 
tween  two  homogeneous  forces  changes  at  a  constant  rate,  the  parameters 
A  and  C  may  be  related  to  those  of  range-dependent  attrition-rate  coeffi¬ 
cients  (see  Section  7). 

3.  PREVIOUS  WORK  ON  DEVELOPMENT  OF  ANALYTIC  SOLUTIONS 

Analytic  solutions  to  Lanchester-type  equations  such  as  (4)  are  impor¬ 
tant,  as  stressed  by  Bonder  (see  pp,  30-31  of  reference  10),  for  developing 
insights  into  the  dynamics  of  combat  by  portraying  the  relation  between 
various  factors  in  the  combat  attrition  process  and  the  surviving  numbers 
of  forces  and  also  for  facilitating  sensitivity  and  other  parametric  analyses 
(see  reference  9).  Additionally,  analytic  results  for  such  variable-coeffi¬ 
cient  formulations  have  also  proven  useful  in  the  study  of  the  optimal  con¬ 
trol  of  deterministic  Lanchester-type  attrition  processes,  see  references  37, 
39,  and  40. 

Oberbeck1281  (also  see  reference  33)  was  apparently  the  first  to  give  re¬ 
sults  via  successive  approximations  for  the  general  case  of  (4).  Subse¬ 
quently,  Bonder  and  Farrell1101  and  Taylor1381  have  given  more  explicit  re¬ 
sults,  which  take  advantage  of  the  structure  of  (4),  and  have  applied  these 
general  results  to  specific  functional  forms  of  attrition-rate  coefficients. 

Koopman  (see  pp.  65-67  of  reference  26)  apparently  first  observed  the 
important  result  (see  Note  4)  that  for  a  constant  ratio  of  attrition-rate  co¬ 
efficients  the  time  solution  to  (4)  takes  a  form  no  more  complicated  than 
that  for  constant  coefficients.  When 

a(t)  =  kah(t),  and  b(t)=kbh(t),  (7) 

where  h(t)  denotes  the  common  time-varying  factor  in  the  two  attrition- 
rate  coefficients,  one  obtains 

x(t)  =  zocosh^(f)  —  yo(ka/kt,ysAt\Yup{t),  (8) 

where  \p{t)  =  (kah^fo*  h(s)  ds.  This  result  is  the  same  as  (3)  except  for  a 
transformation  of  the  time  scale  of  battle.  Subsequently,  Isaacs, 1211  Far¬ 
rell,111’1  and  Taylor  (see  Appendix  D  of  references  34  and  35)  have  inad¬ 
vertently  rediscovered  this  result  in  different  modelling  contexts  (see  ref¬ 
erence  36).  Solutions  for  special  cases  of  (7)  had  been  given  earlier  by 
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Bonder.151  Koopman  also  observed  that  we  again  obtain  a  ‘square- 
law’  relation  ka[x0i  —  x2(t)]  =  kb[y0i—y2(t)}. 

Results  are  considerably  more  complex  for  a  nonconstant  ratio  of  attri¬ 
tion-rate  coefficients.  Bonder  and  Farrell1101  and  Taylor1381  have  given  in¬ 
finite  series  solutions  for  the  linear  attrition-rate  coefficients  with  non¬ 
constant  ratio  (6)  for  the  special  case  of  C  =  0  (see  Note  5).  Additionally, 
Bonder  and  Farrell  (see  pp.  296-367  of  reference  10)  used  analogue  com¬ 
puter  methods  to  parametrically  study  the  dynamics  of  this  combat  situa¬ 
tion.  Taylor1381  has  given  results  for  the  power  attrition-rate  coefficients 
(5),  again  for  C  =  0. 

Thus,  all  the  infinite  series  solutions  noted  above  for  a  nonconstant  ratio 
of  attrition-rate  coefficients  apply  only  to  the  restrictive  case  in  which  C  =  0. 
If  w7e  use  (5)  with  m,n> 0  or  (6),  then  we  are  limited  to  modelling  battles 
in  which  the  fire  effectiveness  of  at  least  one  side  is  initially  zero.  [Within 
the  modelling  context  of  a  constant-speed  attack  (see  Section  7),  the  battle 
must  begin  at  the  minimum  of  the  maximum  effective  ranges  of  the  two 
systems.]  The  purpose  of  this  paper  is  therefore  to  remove  the  requirement 
that  C  =  0  and  to  present  solutions  in  terms  of  new  special  Lanchester  func¬ 
tions  that  can  be  tabulated.  This  is  done  for  power  attrition-rate  coeffi¬ 
cients  with  (in  the  context  of  a  constant-speed  attack)  equal  effective  ranges 
of  the  combat  systems  and  linear  attrition-rate  coefficients  corresponding 
to  opposing  weapons  with  different  effective  ranges.  Additionally,  we  de¬ 
termine  what  functional  forms  for  attrition-rate  coefficients  lead  to  simpli¬ 
fication  of  analytic  results. 

4.  A  MATHEMATICAL  THEORY  FOR  THE  SOLUTION  OF  LANCHESTER-TYPE 
EQUATIONS  OF  MODERN  WARFARE 

In  this  section  w7e  develop  a  generalization  of  the  form  of  the  constant- 
coefficient  solution  (3)  to  the  general  case  of  variable  attrition-rate  coeffi¬ 
cients  (see  Note  6).  The  Lanchester-type  equations  (4)  yield  the  X  force 
level  equation 

d2x/dt2~  { [l/a(t)]da/dt}dx/dt—  a(t)b(t)x  =  0,  (9) 

with  initial  conditions  x(t  =  0)  =  x0,  and  \[\/a{t)]dx/dt\t=o= —yn-  The 
solution  of  (9)  is  given  by 

x(t)=CiX1(t)+CiX2(t),  (10) 

where  \xi(t),  %(£)}  denotes  a  fundamental  system  of  solutions  to  (9)  (see 
p.  119  of  reference  20),  and  Ci  and  C2  are  constants  determined  by  the  ini¬ 
tial  conditions. 

We  now  consider  how7  to  choose  xRf)  and  x2{t)  so  that  (10)  with  Ci  and 
C 2  evaluated  will  be  a  generalization  of  (3).  We  recall  the  w7ell  known 
properties  of  the  hyperbolic  functions  that  appear  in  (3) — that,  for  exam- 
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pie,  d[cosh(a6)3f]/df  =  a(6/a)isinh(a6)If,  where  a  and  6  denote  constants. 
We  choose  Xi (t)  and  x2(t)  so  that 

dxi/dt  =  ka(t)y2,  and  dXi/dt  =  ka(t)yi,  (11) 

where  without  further  specification  k  may  be  any  nonzero  constant.  Simi¬ 
larly,  we  choose  and  y2{t)  to  satisfy 

dy\/dt={\/k)b{t)xi,  and  dy2/dt  =  (l/k)b(t)xi.  (12) 

One  easily  verifies  that  if  X\ (f)  and  x2(t)  satisfy  (11)  and  (12),  then  they 
will  also  satisfy  the  differential  equation  (9).  Next,  we  consider  P(t)  = 
%i(t)yi(t)  —  Xi(t)yt(t).  We  observe  that  dP/dt  is  defined  and  equal  to  zero 
Vfe(f0,  +  °o )  where,  assuming  that  a(t)  and  6(f)  are  both  not  constant, 

fo  =  max(f0A,  f0r).  (13) 

Ux  denotes  the  largest  finite  singularity  (see  p.  160  of  reference  20)  on  the 
f-axis  for  the  X  force-level  equation,  and  similarly  for  tY  (see  Note  7). 
(We  take  to*  to  be  —  oo  if  (9)  has  no  finite  singularity.)  Thus,  for  te 
(to,  +  <» )  we  have  P(t)  =  C s,  where  CN  denotes  a  normalization  constant 
that  we  may  take  to  be  equal  to  one  without  loss  of  generality.  Hence, 
we  can  always  achieve 

xi(t)yi(t)-x2(t)y2(t)  =  l,  (14) 

Vfc(<0,  +  <»  ). 

The  general  Lanchester  functions  xx (t),x2(t),  Vi{t),  and  y2(t)  may  be  con¬ 
structed  by  either  successive  approximations  (see  reference  38)  or  infinite 
series  methods  (see  references  20,  22,  and  27).  Without  further  specifica¬ 
tions,  however,  they  are  not  uniquely  determined.  Moreover,  if  we  can  find 
four  such  functions  satisfying  (11),  (12),  and  (14),  then,  for  example,  x\ (f) 
and  22(f)  are  a  fundamental  system  of  solutions  to  (9). 

We  now  show  that  to  obtain  generalizations  of  the  hyperbolic  functions 
we  must  choose 

k=(kb/ka)K  (15) 

If  a(t)  and  6(f)  are  both  bounded  at  t0,  then  a(f=f0)6(f  =  f0)  =0  and  dP/dt 
=  0  Vfe[fo,  +«).  In  this  case  we  can  achieve  the  normalization  (14) 
Vfe[f0,  4-  °° )  and  uniquely  determine  the  general  Lanchester  functions  by 
choosing 

Xi(t=t0)  =yi(t=t0)  —  1,  and  x2(t  =  f0)  =  2/2  (f  =  to)  =  0.  (16) 

When  (7)  holds  with  h(t  =  0)  =0  and  0<h(t)  <  +  =0  for  <>0,  we  may  use 
Theorem  2  to  solve  (11),  (12),  and  (16).  We  thus  obtain 

21(f)  =  yi(t)  =  cosh0(<),  and  22(f)  =2/2(0  =  sinh0(f),  (17) 

where  6(t)  =  (kakb)i JV  h(s)  ds.  Consequently,  (11)  and  (17)  yield  (15). 


50 


James  G.  Taylor  and  Gerald  G.  Brown 


Next,  we  show  that  if  we  choose  X\ (t)  and  x2(t)  as  above,  then  they  are 
linearly  independent.  Considering  the  Wronskian  of  xh  x2,  denoted  as 
W(x i,  x2),  we  see  that 


W(x i,  x2)  = 


Xi 

Xl 


Xi 

Xi 


a(t)(kb/ka)\ 


Thus,  W(x i,  x2)  9^0  Vte(t0,  +  “  )  (see  Note  8),  and  the  linear  independence 
is  proven. 

We  finally  discuss  how  to  uniquely  specify  Xi(t),  x2(t),  yi(t),  and  y2{t) 
so  that  they  may  be  constructed  by,  for  example,  infinite  series  methods. 
If  Xi(t),  x2(t),  yi(t),  and  y2(t)  are  bounded  at  t  =  t0,  then  they  are,  at  least 
in  a  one-sided  sense,  continuous  at  t  =  t0;  and  we  impose  (16).  If  X\ (t), 
x2 (t),  yi(t),  and  y2(t)  are  not  all  bounded  at  t  =  t0,  we  must  do  other¬ 
wise.  For  the  Fuchsian  equation  (see  p.  370  of  reference  20  and  pp.  108- 


TABLE  I 

Properties  of  the  General  Lanchester  Functions  Xi,  Xi,  yi,  yi 

1 
2 

3 

4 

5 


dxi/dt  =  (WL>)!o(02/2 
dx2/dt  =  (h/ka^ai^yi 
dyi/dt=  (kjkb^bi  t)x% 
dy2/dt=  (ka/kb)ib(t)xi 
x\{t)yi{t)-x2{t)y2{t)  =  l  V< 


127  of  reference  27)  that  arises  from  the  power  attrition-rate  coefficients 
(5)  when  m+n+2  =  0,  we  take  Xi(t  =  tQ)  to  be  unbounded  [so  that  y2(t  =  to) 
is  also].  Then  the  conditions  Zi(f  =  0)  =  l,  (11),  (12),  and  (14)  uniquely 
determine  the  general  Lanchester  functions. 

The  properties  of  the  general  Lanchester  functions  are  summarized  in 
Table  I.  In  general,  the  functions  x-i(t)  and  yi(t)  are  like  the  hyperbolic 
cosine  (with  the  appropriate  argument),  while  x2{t)  and  y2(t)  are  like  the 
hyperbolic  sine.  Using  the  specifications  given  in  the  previous  paragraph, 
we  may  construct  the  general  Lanchester  functions  by  either  successive 
approximations  (see  reference  38)  or  infinite  series  methods  (see  references 
20,  22,  27). 

The  constants  C i  and  C2  in  (10)  are  determined  by  the  initial  conditions 
for  (9),  and  thus  by  (11)  and  (14)  it  follows  that 

x(t)  =  x0{yi(t=0)x1(t)-y2(t=0)x2(t)} 

i  ( lo) 

—  yo(k„/kby{xi(t  =  0)x2(t)  —  x2(t  =  0)xi(t)  j, 

where  a(t)  =  kag(t )  and  b(t)  =kbh(t).  Here  g{t)  and  h(t)  denote  time- 
varying  factors  of  a(t )  and  b(t)  such  that  a(t)/b(t)  =ka/kb  for  g(t )  =h(t). 
Expression  (18)  is  the  generalization  of  the  well  known  constant-coefficient 
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form  (3)  of  the  X  force  level,  x(t),  to  the  general  case  of  variable  Lanchester 
attrition-rate  coefficients.  Equation  (14)  plays  a  key  role  in  verifying 
that  (18)  does  indeed  satisfy  the  initial  conditions  for  (9). 

Recalling  the  well  known  constant  ratio  of  attrition-rate  coefficient  re¬ 
sults  (8),  we  know  that  for  t0^0  (18)  is  at  least  sometimes  capable  of  fur¬ 
ther  simplification.  However,  it  is  indeed  remarkable  that  (18)  can  be 
simplified  only  when  (7)  holds. 

Theorem  1.  For  to  9* 0,  one  can  further  simplify  (18)  if  and  only  if  (7) 
holds  ( constant  ratio  of  attrition-rate  coefficients). 

Proof.  From  (11)  and  (12)  we  have,  for  example,  that  dx\/dy2  = 
k1{a(t)/b{t)}yi/xi.  Thus  a  relation  that  is  independent  of  t  exists  between 
X\  and  y2  if  and  only  if  a (t )/b (f)  =  constant.  Moreover,  when  such  a  rela¬ 
tion  does  exist,  we  have  by  (7)  that 

[zi(0]2-h(0]2=l.  (19) 

Besides  being  a  classic  Lanchester  ‘square  law,’  equation  (19)  is  a  neces¬ 
sary  and  sufficient  condition,  for  example,  for  xi (t)  to  possess  a  so-called 
algebraic  addition  theorem  (see  Note  9).  Unless  the  general  Lanchester 
functions  possess  such  algebraic  addition  theorems,  for  t09^0  there  is  no 
further  simplification  to  (18). 

Remark  1:  When  (7)  does  hold  with  l0< 0,  we  have  (17)  with  B(t)  = 
(kah)ifo‘h(s)  ds+(kakb)if°tl)h(s)  ds.  The  complementary  X-  and  7-func¬ 
tions  are  identical  only  in  this  case.  Moreover,  (14)  [or,  equivalently, 
(19)]  is  then  a  well  known  property  of  the  hyperbolic  functions.  Using 
the  algebraic  addition  theorems  for  the  hyperbolic  functions  [e.g.,  coshit 
coshr  —  sinhwsinhr  =  cosh(u  —  v)],  we  may  simplify  (18)  to  (8).  We  also 
note  that  (8)  may  be  written  as 

x(t)  =  x0cosh([a(t)b(t)]H)  —  y0[a(0/^(0]isinh([a(f)6(f)]i<),  (20) 

where  [a(f)b(f)]J  denotes  an  average,  i.e.,  [a(t)b(t)Y  =  (1  /t)  JY[a(s)6(s)]’  ds. 
Thus,  we  see  that  in  the  case  of  a  constant  ratio  of  attrition-rate  coefficients, 
the  X  force  level,  x(t),  is  given  by  an  expression  formally  equivalent  to 
that  for  the  constant-coefficient  case  with  averages  being  used. 

Remark  2\  In  cases  in  which  a  square-law  relation  like  (19)  does  not  hold, 
there  is  no  such  convenient  reduction  of  (18)  with  0  to  a  simpler  form 
like  (8)  via  an  algebraic  addition  theorem.  Thus,  in  general,  the  X  force 
level,  x(t),  does  not  take  a  simple  form  except  when  f0  =  0  or  a(t)/b(t)  is  a 
constant. 

Remark  3:  When  to  =  0,  then  (18)  again  takes  a  particularly  simple  form: 

x(t)  =  XoX!(t)  —yo(ka/kb)ix2(t) .  (21) 

Our  previous  results  in  reference  38  were  all  of  the  form  (21). 

In  reference  38  (see  also  references  35  and  36)  we  showed  that  when  the 


52 


James  G.  Taylor  and  Gerald  G.  Brown 


ratio  of  attrition-rate  coefficients  is  constant,  we  can  transform  the  X  force- 
level  equation  into  one  with  constant  coefficients  by  a  transformation  of  the 
independent  variable  t.  As  we  have  seen,  this  case  leads  to  particularly 
convenient  results.  In  this  respect,  a  useful  theorem  (see  also  Appendix 
D  of  reference  34)  is 

Theorem  2.  A  necessary  and  sufficient  condition  to  be  able  to  transform  the 
X  force-level  equation  (9)  by  a  transformation  of  the  independent  variable  t 
into  a  linear  second-order  ordinary  differential  equation  with  constant  coeffi¬ 
cients  is  that  {[l/a(t)}da/dt—[l/b(t)]db/dt\/[a(t)b(t)Y  is  a  constant. 

Proof.  The  theorem  follows  immediately  from  a  result  given  on  pp.  73- 
74  of  reference  4:  the  equation  d2x/dt2-\-p(t)dx/dt-\-q{t)x  =  0  can  be  trans¬ 
formed  to  a  linear  equation  with  constant  coefficients  by  a  transformation 
of  the  independent  variable  alone  if  and  only  if  the  expression  {{dq/dt+ 
2pq)/q312}  is  found  to  be  constant. 

Moreover,  when  the  X  force-level  equation  is  so  transformable,  the  de¬ 
sired  substitution  is  given  by  u  =  K  j‘  [a(s)b(s)Y  ds,  where  /*•  •  ■  ds  denotes 
an  indefinite  integral  and  K  is  an  arbitrary  constant  conveniently  chosen. 

For  example,  for  the  power  attrition -rate  coefficients  (5),  Theorem  2  tells 
us  that  we  can  transform  (9)  into  an  equation  with  constant  coefficients 
only  when  m  =  nor  m+n-Y 2  =  0. 

5.  POWER  ATTRITION-RATE  COEFFICIENTS 

In  this  case  the  attrition-rate  coefficients  are  given  by  (5).  There  are 
two  cases  to  be  considered,  depending  on  whether  (9)  with  (5)  can  be  trans¬ 
formed  into  a  constant-coefficient  equation  (see  Theorem  2  above):  (1) 
m+n- 1-2^0,  and  (2)  m  +  n-\- 2  =  0. 

Case  1.  m+n+ 2^0 

In  order  for  solution  methods  (either  successive  approximations  or  in¬ 
finite  series  methods)  to  be  applicable,  we  must  further  impose  the  following 
restrictions:  (a)  for  C  =  0,  we  must  have  rn >  —  1  and  n>  —  1,  while  (b) 
for  C > 0,  we  must  have  rn  +  n  +  2^0.  Only  the  former  case  was  considered 
in  reference  38.  The  theory  presented  in  the  previous  section  is  essential 
for  extending  these  results  to  the  case  in  which  C  >  0. 

From  results  given  in  reference  38,  it  follows  that  a  fundamental  system 
of  solutions  to  (9)  with  attrition -rate  coefficients  (5)  is  given  by 

^i(0  =  r(?)t(A:oA:4)V2s]3,(f+C')(m+1),2/_p[(A:A)i(f+0)7s],  (22) 

and 

x*(t)  =  r(p)[(fcoUV2s]s(f+C)<m+1),7P[(fcafc6)7f+0)7s],  (23) 

where  Ip  denotes  the  modified  Bessel  function  of  the  first  kind  of  order  p, 
p=(m+l)/(m+n+2),  p+q=  1,  and  2s=m+n+2,  with  similar  results 
holding  for  \yi(t),  yi(t)}.  The  above  Lanchester  functions,  of  course,  sat¬ 
isfy  all  the  properties  given  in  Table  I.  The  X  force  level,  x(t),  then  is 
given  by  (18). 
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To  emphasize  the  dependence  on  the  parameter  C  it  is  sometimes  con¬ 
venient  to  write,  for  example,  xi(t)  =  Xi(t;  C ).  For  a(t)  and  b(t)  both  not 
constant,  the  only  finite  singularity  of  the  X  and  Y  force-level  equations 
with  the  power  attrition -rate  coefficients  (5)  is  t=  —C  so  that  t0=—C. 
We  note  that  then  xi(t  =  0;  C  =  0)  =  y1(t  =  0;  (7  =  0)  =  1  and  Xz{t  =  0;  (7  =  0) 
=  yi{t  =  0;  (7  =  0)  =0;  hence  (18)  with  these  particular  fundamental  sys¬ 
tems  reduces  to  (20)  when  (7  =  0.  By  Theorem  1  no  such  simplification  of 
(18)  occurs  for  C>0  unless  m  =  n. 

The  above  representations  (22)  and  (23)  of  fundamental  systems  of  solu¬ 
tions  to  the  X  and  Y  force-level  equations  are  not  particularly  convenient. 
For  example,  for  m,  n>  —  1  wehaveOCp,  q<  1,  and  p+q  =  l,and  tabulations 
only  exist  of  the  modified  Bessel  function  of  the  first  kind  of  fractional  order 
v  for  a  restrictive  set  of  values  [i.e.,  v=  ±}/i,  (see  reference 

1)].  This  Bessel  function  reduces  to  other  tabulated  forms  (i.e.,  hyberbolic 
functions)  for  Therefore,  it  is  more  useful  to  express  x\ (t), 

xi (t),  yi(t),  and  y^t)  in  the  form  of  infinite  series  and  to  consider  the  result¬ 
ing  transcendental  functions  as  entities  in  their  own  right.  Thus,  it  is  con¬ 
venient  to  define  the  following  Lanchester  functions  for  power  attrition- 
rate  coefficients  (referred  to  as  power  Lanchester  functions) . 

=  r (q)  Zt-o  [(kah)*/m+n+2)Tktklm+n+i)/ 

[k\T(k  +  q)], 

»».•«)  =  r(p)  E£?[(W t/(m+n+2)fk+1tUm+n+2)+m+l/ 

[klT(k+l+p)], 

Um,n(t)  =  r(p)  [(kakb)i/(7n+n+2)fktHm+n+2)/ 

[k'.T(k+p)}, 

VmAt)  =  r(g)  Lto”  [(kMi/(m+n+2)]u+1tHm+n+2)+n+1/ 

[k\V(k+l+q)}. 


(24) 


(25) 


(26) 


(27) 


The  power  Lanchester  functions  um,n,  vm,n,  Um,n,  and  Vm,n  have  the  prop¬ 
erties  shown  in  Table  II.  For  reasons  to  be  explained  below  we  will  refer 
to  um,n(t)  and  Um,n(t)  and  similarly  vm,n(t)  and  Vm,n(t)  as  complementary 
power  Lanchester  functions.  Property  5  of  Table  II  follows  from  proper¬ 
ties  6  and  7  for  m,  n>  —  1.  We  show  below  that  it  holds  for  all  values  of 
m  and  n.  Thus,  for  (7^0,  with  m  and  n  subject  to  the  restrictions  given  in 
the  paragraph  immediately  preceding  equation  (22),  the  solution  (18)  to 
(9)  with  the  power  attrition -rate  coefficients  (5)  may  then  be  written  in 
terms  of  the  power  Lanchester  functions  as 


X(t)—Xo{  L  m,n(C)um,n(t-{-C)  —  V m,n(C)vm,n(t-\-C)  } 

Voika/kb)  {'amfn{C)V7n,n{t~\~C')  Wmin(t-\-C')Vmin(<C')  }  , 


(28) 


and  similarly  for  y(t).  From  (28)  we  see  that  the  methods  of  reference  38 
would  become  hopelessly  bogged  down  in  details  for  C>0. 
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TABLE  II 


Properties  of  the  Power  Lanchester  Functions  um,n,  vm,n 

j  Um,ny 

AND  Vm,n 

1. 

dUm  ,n / dt  =  {kb/ ka)^[katm\  7m  ,n  (0 

2. 

dvm,n/dt=  {kb/  ka)^[katm]  U  m  ,nif) 

3. 

dUm,n/dt=  (ka/kb)^[kbtn]Vm,n{i) 

4. 

dVm,n/dt  —  (ka/ kb)^[kbln]um  ,n  (0 

5. 

Pm  ,n(0 Cm  ... (0  Pm,n(0  F m,nd  =  1  V i 

6. 

Pm,n(<=0)  =  f7m,»(2  =  0)  =  l 

7. 

Wm,n(«==0)  =  7m,n(<  =  0)=0 

8. 

um,m(t)  =  Um,m(t)=  eosh[  (fcaA:i,)hm+1/(w+ 1)] 

9. 

Pm,  m(i)  =  7m, rn  (0  =sinh[  ( kah)Hm+ V  (m+ 1)1 

Note 

:  Properties  6  and  7  hold  only  for  m,»>  —1. 

For  computational  reasons  it  is  convenient  to  introduce  the  auxiliary 

power 

Lanchester  functions  [also  referred  to  as  the  Lanchester-Clifford- 

Schlafli  functions,  or  LCS  functions  (see  Note  10)],  which  are  defined  for 

p^O,  - 

-1,  -2,  —3,  •  •  • 

F,{x)=YXZ 0  (x/2)n/[kiuk-l  (j  +  p)], 

(29) 

and 

Gy(x)=Zt-o  (x/2)u+1/[klri-k0  (j  +  p)], 

(30) 

where  we  have  adopted  the  convention  that  11;=* /,-  =  1  for  N<M.  The 
LCS  functions  possess  the  properties  shown  in  Table  III.  Considering 
Q(x)  =  Fp(x)Fq(x)—Gp(x)Gq(x),  where  p+q=  1  and  p,  q^O,  —1,  —  2,  •  •  •, 
one  may  easily  prove  property  3  of  Table  III  by  observing  that  dQ/ dx= OVx 
and  using  properties  4  and  5. 

The  power  Lanchester  functions  um,n,  vm,n,  Um,n,  and  Vm,„  may  be  ex¬ 
pressed  in  terms  of  the  LCS  functions  as  follows: 

V»W=W(i)],  (31) 

TABLE  III 

Properties  or  the  LCS  Functions  F„(x)  and  G,(x). 


1.  dF„/dx=G,(x ) 

2.  dG„/dx  =  F  ,{x)—\(2v—l)/x\G,{x) 

3.  F„(x)Fq{x)-GP(x)Gq(x)  =  1  Vx 

where  p+2  =  l  and  p,q^0,  —1,  —2,  —3,  ••• 

4.  F,(x=0)  =  l 

5.  Gy(x=Q)  =0 

6.  limI^o  {G,(x)/x}  =  l/(2r) 

7.  dFy/dx{x  =  0)=0 

8.  dG,/dx{x= 0)  =  l/(2v) 

9.  Fi/s(x)  =  coshz 

10.  Gi/8(x)=sinhx 
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vm,n(t)=tim~n)l2Gp[Ht)],  (32) 

Um,n{t)  =  FPm)],  (33) 

=  (34) 


where  \p{t)  =  {[katm][k}>tn\),!t/[{m+n  +  2)/2].  Noting  that  um,n{t)U m,n(t) — 
vm,n(t)VmAt)  =Fp[ip(t)]F<l[ip(t)]  —  Gp[ip(t)]Gv[t(t)],  we  see  that  property  5 
of  Table  II  holds  for  all  values  of  m  and  n  and  follows  from  property  3  of 
Table  III.  Also,  um,n(t )  and  Um<n(t)  have  been  called  complementary 
power  Lanchester  functions  because  of  the  above  relations,  the  fact  that 
p+q  =  1,  and  property  3  of  Table  III.  Consequently,  (28)  may  now  be 
written  as 

x(t)  =  xo\FAe(0)]FMt)}-(i+t/cym-n)liGMO)}GP[em 

-yola(t  =  0)/b(t  =  0)\H(l+t/Cym-n)l2FMO)}GPW)]  (35) 

-Gp[0(O)]F9[0(O]}, 

where  p=  (m+l)/(m+n+2),  q=l—p,  and  9(t)  =  [a(i!)&(i!)]s(<+C)/ 
[{m+n+ 2)/2],  Introduction  of  the  LCS  functions  thus  sheds  light  on  the 
parametric  dependence  of  solutions:  there  are  two  exponent  parameters,  p 
and  (m  —  n)/ 2,  and  an  ‘intensity’  parameter,  9(t),  as  well  as  the  relative 
effectiveness  parameter  ( ka/kb)K  We  observe  that  average  ‘intensity’ 
[a(f)6(i)]J  is  related  to  6(t)  =  9(t;  C)  by  [a(t)b(t)]H=  6(t\  C)  —  6(t  =  0;  C). 
Without  an  algebraic  addition  theorem,  however,  this  does  not  lead  to 
simplification  of  (35).  We  do  find,  though,  that  for  C  =  0  with  m,  n>  —1 
we  have  [a(t)b(t)]h  =  9(t;  C  =  0),  and  (35)  may  be  written  like  (21)  as 

x(t)  =  x0Fq{[a{t)b{t)]ht)  -y0[a(t) /b{t)]hGp{[a{t)b{t)]H).  (36) 

It  may  be  shown  that  lim(_o  [a(0/f>(<)]*C!,([a(<)6(f)]if)  =0. 

Considering  (35)  and  the  above,  we  see  that  there  is  an  interesting  way 
(with  some  similarities  with  the  constant  coefficient  situation)  of  thinking 
about  combat  between  two  homogeneous  forces  described  by  (4)  and  (5) : 
the  course  of  combat  is  determined  by  average  combat  ‘intensity’  and  howr 
it  changes  over  time  as  well  as  the  initial  force  levels,  xo  and  yo,  and  [a(t  = 
0)/b(t  =  0)}\ 

Case  2.  m+n+2  =  0  and  C>0 

It  should  first  be  noted  that  the  solution  obtained  for  Case  1  becomes  in¬ 
determinate  [consider,  for  example,  what  happens  to  um,n(t)  as  defined  by 
(24)].  To  solve  the  Fuchsian  equation  (see  pp.  108-127  of  reference  27) 
which  arises  from  (9)  with  (5)  wffien  w+n+2  =  0,  we  make  the  substitu¬ 
tion  «  =  ln(f+C)  to  transform  (9)  into  an  equation  with  constant  coeffi¬ 
cients.  Solving  this  equation  and  transforming  back  to  the  original  inde¬ 
pendent  variable  t,  we  find  that  both  members  of  the  resulting  fundamental 


56 


James  G.  Taylor  and  Gerald  G.  Brown 


system  of  solutions  are  not  bounded  at  t  — 10=  —C  (see  pp.  148-150  of  refer¬ 
ence  29),  and  we  consequently  cannot  impose  (16).  To  determine  a 
fundamental  system  uniquely,  we  take  Xi(t=t0)  to  be  unbounded  and 
specify  the  conditions  £i(f  =  0)  =  l,  (11),  (12),  and  (14)  (see  Section  4). 
The  corresponding  fundamental  system  (with  the  properties  given  in  Table 
I)  for  (18)  is  given  by 

*i  (t)  =  ([t+C)/C)a-,  and  x2(t)={C-ln+1)/2y}(a-d-)H[t+C]/Cr+, 

and  similarly  for  {yi(t),  where  y  —  {kakb+[(ni  +  l)/2f\ \  a+  =  7  + 

(m  + 1)/2=  —  /3-,  and  a~=  — 7+(m4T)/2  =  —  /3+.  In  this  case  (18)  be¬ 
comes 

x«)  =  Xo![l/(27)][^+([<+C]/C)“+-d-([<+C]/C)“-} 
-^(fca/^)i{[(U6)^m+1/(27)][([i+C]/C)“+-([«+C]/C)a-]!, 
which  may  also  be  written  in  the  more  convenient  form 

m  =  at+C}/C(m+1)i2{x0[cosh[y\n([t+C]/C)] 

-[(OT  +  l)/(27)]sinh(7[f+C]/C)]  (38) 

- yo(ka/kb)  \{Cm+\kakb)  V7)sinh{7ln(  [t+C]/C) }  ] } . 

Although  the  above  solutions  appear  complicated,  they  are  readily  evalu¬ 
ated  with  the  help  of  a  portable  electronic  calculator  such  as  is  commercially 
available  today.  In  fact,  such  a  calculator  can  even  be  ‘programmed’  to 
facilitate  parametric  analyses. 

6.  WEAPON  SYSTEMS  WITH  DIFFERENT  EFFECTIVE  RANGES: 

LINEAR  ATTRITION-RATE  COEFFICIENTS 

Another  interesting  situation  in  which  to  apply  our  general  theory  is  that 
of  combat  between  two  homogeneous  forces  that  use  weapons  with  different 
effective  ranges.  We  consider  the  example  previously  examined  by  Bonder 
and  Farrell1101  and  Taylor1381  of  a  constant-speed  attack  of  a  mobile 
force  against  a  static  defense.  The  weapon  systems  of  the  two  sides  have 
different  effective  ranges,  and  the  lethality  of  each  side’s  fire  depends  line¬ 
arly  upon  range.  We  assume  that  the  opening  range  of  battle,  Ro,  satisfies 
Ro^ minimum  (Ra,  Rp),  where  Ra  denotes  the  maximum  effective  range  of 
the  F  system.  The  Lanchester  attrition-rate  coefficients  for  such  a  battle 
may  be  written  as  (6).  The  parameter  C  is  related  to  the  opening  range  of 
battle  in  comparison  with  minimum  (Ra,  Rp),  whereas  A  reflects  the  dif¬ 
ference  in  maximum  effective  ranges  [see  (51)  in  Section  7].  Figure  1 
shows  the  range  dependencies  of  these  Lanchester  attrition-rate  coefficients 
and  the  opening  range  of  battle. 

From  results  given  in  reference  38,  it  follows  that  a  fundamental  system 
of  solutions  to  (9)  (which  has  the  properties  shown  in  Table  I)  with  the 
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attrition -rate  coefficients  (6)  (and  also  one  for  the  Y  force-level  equation) 
is  given  by 

xl(t)=f(t+C),  x2(t)  =  g(t+C), 

yi(t)=F(t+C),  and  y^t)  =G(t+C), 

f(t)  =  i[(kakb)i/2fn/(2n) !}  ECo"  BnkAktin-k,  (40) 

g(t)=Zi: o  U(fcA)V2]2”+1/(2n+l) '}  Eto”  CnkAkt4n+2~k,  (41) 


Fig.  1.  Linear  attrition-rate  coefficients  for  weapon  systems  with  different  ef' 
fective  ranges.  [Notes:  1.  The  maximum  effective  ranges  of  the  two  weapon  sys¬ 
tems  are  denoted  as  Ra  and  Rp.  2.  The  opening  range  of  battle  is  denoted  as  Ro 
and  (as  shown)  #0<minimum  (Ra,  Rp).] 


m  =  E:=“  {[( W i/2f7(2n) !}  Et?  DnkAktin~k,  (42) 

G(t)  =  E"=?  \{(kakb)i/2fn+1/(2n+ 1) !}  E?=o  EnkAkt4n+2~k.  (43) 


The  coefficients  Bnk  and  Cnk  are  given  in  reference  38,  while  the  coefficients 
Dk  and  Ek  are  given  by  D0°  =  1  and  for  n> 0 


1, 

[4«(4n  —  2)/(4n  —  k)]{[Dkn-if  (4n  —  k  —  2)} 

+  [Dt-\/(4:n-k-l)]}, 


for  k  —  0, 

for  lf£/c^n— 1, 
for  k  =  n; 


{(y3)[(in-2)/(3n-l)]DnnZ\, 
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and  E0°=  1,  Ef  =  2,  and  for  n> 0 

fl,  for  k  =  0, 

„  I  [4n(4n+2)/ (4n  — fc+2)]{[-E't_i/(4n  — /c)] 

E  w  —  < 

+  [-Bb-i/(4/i— fc  +  1)]},  for 

,(M)[(4w+2)/(3n+l)]f?^_x,  for  fc  =  m+l, 

where  we  have  adopted  the  convention  that  0°  =  1 .  We  note  that  f{t  =  0 )  = 
F(<  =  0)  =  1  and  g(t  =  0)  =  G(t  =  0)  =  0.  Motivated  by  their  origin  within 
the  context  of  a  constant-speed  attack  of  mobile  forces  as  first  studied  by 
Bonder/5,7 ,10i  we  will  refer  to  the  above  Lanchester  functions  for  the  linear 
attrition-rate  coefficients  with  nonconstant  ratio  (6)  as  offset  linear  Lan¬ 
chester  functions.  We  use  the  word  ‘offset’  to  denote  a  nonzero  difference 
between  the  maximum  effective  ranges  of  the  two  systems  (consider  Fig. 
1). 

To  emphasize  the  dependence  on  the  parameter  4,  it  is  sometimes  con¬ 
venient  to  write,  for  example,  f(t)=f(t;  A).  We  note  that  /(i;  .4  =  0)  = 
F(t;  4  =  0)  =  cosh[(fcBfcj)  V/2]  and  g(t;  A  =  0)  =  G(t;  4  =  0)=sinh[(fcofa>)i<2/2]. 
For  4  =  0,  the  fundamental  property  f(t)F(t)— g(t)G(t)  =  l  of  these  offset 
linear  Lanchester  functions  becomes  a  well  known  property  of  the  hyper¬ 
bolic  functions. 

As  was  the  case  for  the  power  attrition-rate  coefficients  of  the  previous 
section,  for  computational  reasons  it  is  convenient  to  introduce  the  follow¬ 
ing  auxiliary  offset  linear  Lanchester  functions  (depending  on  two  parame¬ 
ters)  : 

{\U/(2n)l}Zl-S  Bny,  (44) 

w(\,  y)  =  E:=o”  { X2”+1/ (2n+l) !}  DJ: 0  C«V,  (45) 

H(\,  n)  =  Zl-t  {A27(2 n) !}  Zk-S  Dn V,  (46) 

W(\,  y)=Znn=o  {X2"+1/ (2n+ 1) !}  7V-  (47) 

The  auxiliary  offset  linear  Lanchester  functions  possess  the  following  prop¬ 
erties  : 

1.  h(\,  fi)  —  w(\,  y)W(\,  u)  =  1  VX,  u, 

2.  h(\,  0)  =H(\,  0)  =  coshX, 

3.  io(X,  0)  =  TF(X,  0)  =  sinhX. 

Using  the  auxiliary  offset  linear  Lanchester  functions,  we  may  write  the 
X  force  level  as 

x(t)  =Xo{H[0(O),  8(0)]A[8(0,  ««)]- W[e{Q),  «(0)M«(0,  8(01} 
-Vo(ka/h)Hh[0(0),  8(0)M«(0,  8(0] 

-iw[e(0),  5(O)]A[0(O,  8(0]}) 


(48) 
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where  S(t)=A/(t+C)  and  d(t)  =  (kakby(t+C)2/ 2.  The  expression  S(t) 
reflects  how  much  the  above  Lanchester  functions  deviate  from  the  hyper¬ 
bolic  functions. 


7.  SOME  NUMERICAL  EXAMPLES 

In  this  section  we  examine  three  numerical  examples  that  illustrate  possi¬ 
ble  use  of  some  of  our  new  results.  These  examples  are  extensions  of  the 
ones  given  previously  in  reference  38  by  having  the  opening  range  of  battle 
be  less  than  the  minimum  of  the  two  maximum  effective  ranges  for  the  two 
weapon  systems.  They  are  motivated  by  the  work  of  Bonder16,7,101  on  the 
value  of  range  capabilities  and  mobility  for  weapon  systems  in  combat  de¬ 
scribed  by  Lanchester-type  equations  of  modern  warfare.  The  modelling 
context  of  these  examples  is  that  of  weapon  systems  with  (a)  different 
range  dependencies  of  lethality  of  each  side’s  fire,  but  the  same  maximum 
effective  range,  and  (b)  linear  attrition-rate  coefficients,  but  different  ef¬ 
fective  ranges. 

As  in  reference  38  we  consider  a  constant-speed  attack  on  a  static  defen¬ 
sive  position  with  the  combat  dynamics  described  by 

dx/dt=  —a(r)y=  -a0(l -r/Ra)my, 

(49) 

dy/dt—  ~p(r)x=  ~ft0(l  —  r/Ri})nx, 

where  Ra  and  R$  denote  the  maximum  effective  ranges  of  the  Y  and  X 
weapon  systems,  respectively  (i.e„  a(r)  =0  for  r>Ra).  The  reader  should 
consider  (49)  to  be  (4)  with  the  attrition-rate  coefficients  (5)  or  (6)  (as  is 
appropriate)  expressed  in  terms  of  range  (i.e.,  force  separation).  Range  is 
related  to  time  by  r(t)  —R0—vt,  where  R0  denotes  the  opening  range  of  bat¬ 
tle.  Several  range  dependencies  for  an  attrition-rate  coefficient  are  shown 
in  Fig.  2,  and  an  opening  range  less  than  the  weapon  system’s  maximum 
effective  range  is  indicated.  The  parameters  of  the  attrition-rate  coeffi¬ 
cients  ce[r(t )]  and  /3[r(<)]  in  (48)  are  readily  related  to  those  of  a(t)  and  b(t) 
in  (5)  and  (6).  For  (5)  we  have 

ka  =  ao(v/Ra)m,  kb  =  Po(v/Rp)n,  and  C=(Ra—R0)/v,  (50) 

where  Ra  =  R^Ro,  while  for  (6)  we  have 

ka  =  ctov/Ra,  kb=0ov/Rp,  A  =  (Rs—Ra)/v,  and  C=  (Ra— Ro)/v,  (51) 
where  Rp^Ra^Ro- 

Numerical  results  are  shown  in  Figs.  3,  4,  and  5.  The  force-level  tra¬ 
jectories  have  been  generated  by  a  digital  computer  program  using  the  aux¬ 
iliary  Lanchester  functions  [i.e.,  F,  and  G,  for  (35)  and  h,  w,  H,  and  W  for 
(48)].  These  functions  are  particularly  convenient  for  such  digital  com¬ 
puter  work.  The  computer  routines  were  checked  against  the  numerical 
results  given  previously  in  reference  38.  In  this  work  we  have  taken  the 
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opening  range  of  battle  R0  <  minimum  (Ra,  R^).  Numerical  values  for 
battle  parameters  (except  those  for  R0,  Ra,  and  Rs)  are  the  same  as  those 
used  for  the  examples  considered  in  reference  38. 

For  Figs.  3  and  4  both  weapon  systems  have  the  same  maximum  effective 
range  (i.e.,  Ra  =  Ra).  As  done  for  the  plots  in  reference  38,  we  have  held 
«o  =  a(r  =  0)  and  do  constant  and  have  varied  the  exponents  m  and  n,  which 


range  /i  (meters) 

Fig.  2.  Dependence  of  the  attrition-rate  coefficient  a(r)  on  the  exponent  m  for 
constant  maximum  effective  range  of  the  weapon  system  and  constant  kill  capa¬ 
bility  at  zero  range.  [Notes:  1.  The  maximum  effective  range  of  the  system  is  de¬ 
noted  #a=2000  meters.  2.  a(r=0)  =a0=0.6  X  casualties/(unit  time  X  number 
of  F  units)  denotes  the  Y  force  weapon  system  kill  rate  at  zero  force  separation 
(range).  3.  The  opening  range  of  battle  is  denoted  as  #0  =  1250  meters  and  (as 
shown)  #o  <#«]■ 

control  the  range  dependencies  of  a(r)  and  j3(r).  With  the  exception  of 
#o,  all  the  battle  parameters  for  the  curves  shown  in  Fig.  3  are  the  same  as 
those  used  for  the  corresponding  example  previously  considered  in  reference 
38  (i.e.,  Ra  =  Rp  =  2000  meters).  Consequently,  corresponding  force-level 
trajectories  are  similar,  with  greater  ‘separation’  shown  here  between  curves. 
In  Fig.  3  with  the  opening  range  of  battle  Ro=  1250  meters,  the  curves  cor¬ 
responding  to  the  constant-coefficient  case  (i.e.,  m~n  =  0)  are  exactly  the 
same  (for  the  same  time  intervals)  as  those  shown  in  reference  38  with 
Ro  =  2000  meters.  Other  battle  trajectories  with  m,  n>  0  decay  faster  than 
they  did  in  reference  38  because  the  ‘intensity’  of  combat  is  greater;  i.e., 
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250  500  750  1000  1250 


Fig.  3.  Force-level  trajectories  of  X  and  Y  forces  for  various  combinations  of 
the  exponents  m  and  n  in  the  power  attrition-rate  coefficients  for  Ro  =  1250  meters, 
Ra  =R-j  =2000  meters,  ao=0.06X  casualties/(minutes-  F  unit),  do  =0.6  Y  casual¬ 
ties/ (minutes -X  unit),  e  =  5  mph,  a;0  =  10,  and?/o  =  30.  The  exponent  combinations 
are  denoted  as  m  :n  in  the  figure,  and  the  symbol  X  denotes  the  end  of  a  force-level 
trajectory  due  to  the  annihilation  of  the  enemy  force. 

as  a  function  of  time  the  attrition-rate  coefficients  are  larger  here  than  in 
reference  38  (see  Fig.  2). 

As  noted  in  references  5  and  38,  knowledge  of  the  range  or  time  de¬ 
pendence  of  weapon-system  kill  capability  is  essential  for  forecasting  the 
battle’s  outcome  from  the  initial  trend  of  battle.  For  example,  compare  in 
Fig.  3  the  outcomes  for  curves  denoted  as  1:0,  1:1,  and  1:2.  Clark11'1 
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has  developed  methodology  for  estimating  such  capability  from  the  output 
of  a  high-resolution  Monte  Carlo  combat  simulation.  Figure  4  shows  similar 
battle  curves  for  the  same  parameter  values,  except  that  Ra  =  Rff=1500 
meters.  Observing  that  for  m^l  we  have  a(r;  Ra)<a(r;  Ra)*=>Ra<Ra, 


250  500  750  1000  1250 


Fig.  4.  Force-level  trajectories  of  X  and  Y  forces  for  various  combinations  of 
the  exponents  m  and  n  in  the  power  attrition-rate  coefficients  for  the  same  parameter 
values  chosen  for  Fig.  3  except  that  Ra  =#3  =  1500  meters.  The  symbol  conven¬ 
tions  are  also  the  same  as  in  Fig.  3. 

we  may  consider  that  the  combat  is  less  ‘intense’  for  such  pairs  of  battle 
trajectories  in  Fig.  4  than  for  those  shown  in  Fig.  3. 

We  consider  next  numerical  results  for  the  linear  attrition-rate  coeffi¬ 
cients  with  nonconstant  ratio  (6).  In  Fig.  5  we  show  the  effect  of  increas- 
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ing  the  effective  range  of  the  defender’s  weapons.  (The  X  force  may  be 
considered  to  be  the  defender.)  For  these  computations,  using  (48)  and 
the  auxiliary  offset  linear  Lanchester  functions,  we  have  held  the  opening 
range  of  battle  constant  at  f?0=1250  meters  and  the  maximum  effective 


250  500  750  1000  1250 


Fig.  5.  Force-level  trajectories  of  X  and  Y  forces  for  various  effective  ranges 
Rg  of  the  X  force  weapons  with  linear  attrition-rate  coefficients  for  Ra  =  1500  meters 
and  the  same  values  of  the  other  parameters  listed  in  the  legend  of  Fig.  3.  The 
symbol  X  has  the  same  meaning  as  in  this  figure. 

range  of  the  Y  weapon  constant  at  Ra  =  1500  meters.  As  in  reference  38, 
both  attrition -rate  coefficients  depend  linearly  on  range  [i.e.,  m  =  n  —  l  in 
(49)],  q:o  and  do  have  been  held  constant,  and  Rg  has  been  varied.  As  shown 
in  Fig.  5,  we  quantiatively  see  the  benefit  from  increasing  the  long-range 
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kill  capability  of  the  defender’s  weapons.  (See  pp.  296-367  of  reference  10 
for  extensive  parametric  studies  based  on  analogue  computer  results.  The 
analytic  results  given  here  may  be  used  for  such  investigations. )  It  should 
be  emphasized  that  although  the  results  shown  in  Figs.  3,  4,  and  5  are 
similar  to  those  in  reference  38,  our  previous  analytic  results  were  limited  to 
the  case  in  which  the  opening  range  of  battle  R0  =  minimum  (Ra,  Rs). 

8.  DEVELOPMENT  OF  TABULATIONS 

Tabulations,  such  as  those  given  in  reference  1,  of  the  auxiliary  Lan- 
chester  functions  (i.e.,  the  LCS  functions  and  the  auxiliary  offset  linear 
Lanchester  functions)  would  lead  to  analysts  being  able  to  generate  nu¬ 
merical  solutions  to  the  variable-coefficient  Lanchester-type  equations  con¬ 
sidered  above  with  somewhat  the  same  facility  as  for  the  constant-coeffi¬ 
cient  case.  One  may  develop  tabulations  of  the  LCS  functions  I'\  ( x )  and 
G„(x)  directly  from  their  power  series  representations  (29)  and  (30). 
In  digital  computer  work  we  have  found  it  convenient  to  calculate  F,{x) 
from  Fy(x)=J2k=ofk(x),  where  fk(x)  is  given  recursively  by  fk(x)  = 
(x/2)2fk~i(x)/[k(k—l  +  v)]  for  k=  1,  2,  •  ■  -.with  fo(x )  =  1  and  similarly  for 
Gv(x).  For  m,  n=0,  1,2,3  in  the  power  attrition-rate  coefficients  (5), 
one  has  the  requirement  for  tabulations  of  F„(x)  and  G>(z)  for  v  =  }i,  %, 
H,  H,  H,  H,  H,  %,  Yi,  and  One  may  use  (14)  to  reduce  the  number 
of  series  evaluations  in  such  work.  Tabulations  of  the  auxiliary  offset  linear 
Lanchester  functions  are  readily  developed  from  the  series  (44)  through 
(47)  with,  for  example,  Dnk  computed  recursively.  We  have  followed  the 
above  computational  procedures  in  the  work  reported  in  Section  7. 

The  determination  of  exactly  which  Lanchester  functions  should  be  tabu¬ 
lated  depends  on  which  functional  forms  for  attrition-rate  coefficients  are 
most  important  in  applications,  and  consideration  of  this  is  beyond  the 
scope  of  this  paper.  The  specific  time-dependent,  attrition-rate  coeffi¬ 
cients  we  have  considered  have  been  suggested  by  Bonder’s  work  on  the 
range  dependency  of  Lanchester  attrition -rate  coefficients  (see  pp.  196-200 
of  reference  10). 


9.  SUMMARY 

We  have  presented  a  mathematical  theory  (including  new  general  ana¬ 
lytic  results)  for  the  solution  of  variable-coefficient  Lanchester-type  equa¬ 
tions  of  ‘modern  warfare’  for  combat  between  two  homogeneous  forces. 
These  results  allow  one  to  extend  results  given  in  references  10  and  38  that 
applied  under  rather  restrictive  conditions,  e.g.,  opening  range  of  battle 
equal  to  the  minimum  of  the  two  maximum  effective  ranges  of  the  weapon 
systems.  It  has  been  shown  that  in  general  the  deterministic  time  histories 
of  the  X  and  Y  force  levels  [i.e.,  x(t)  and  y(t)]  may  be  expressed  in  terms 
of  four  complementary  Lanchester  functions  and  not  two  as  in  the  constant- 
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coefficient  case.  The  mathematical  properties  of  these  general  Lanchester 
functions  that  facilitate  analytic  solution  have  been  given.  We  have  paid 
particular  attention  to  the  determination  of  the  conditions  on  attrition-rate 
coefficients  that  allow  relatively  simple  analytic  results.  We  have  applied 
these  general  results  to  two  types  of  attrition-rate  coefficients:  power  at¬ 
trition-rate  coefficients  and  linear  attrition-rate  coefficients  with  noncon¬ 
stant  ratio. 

These  results  may  be  used  to  facilitate  parametric  analysis  of  the  dy¬ 
namic  combat  interactions  between  two  homogeneous  forces  with  time-  or 
range-dependent  weapon  system  capabilities.  Such  models  are  of  particu¬ 
lar  interest  in  light  of  the  work  of  Bonder16'81  and  others  in  references  3,  10, 
30,  and  31  on  the  prediction  of  Lanchester  attrition-rate  coefficients  from 
weapon-system  performance  data  and  the  work  of  Clark1121  on  the  estima¬ 
tion  of  such  time-dependent  coefficients  from  Monte  Carlo  simulation  out¬ 
put.  Additionally,  our  new  theoretical  results  may  be  used  to  solve 
Grubbs  and  Shuford’s  new  probabilistic  formulation  for  Lanchester  combat 
theory  [see  equations  (14)  and  (15)  of  reference  17].  A  further  discus¬ 
sion  of  applications  is  found  in  references  10  and  38. 

NOTES 

1.  This  is  due  to  the  complexity  of  obtaining  analytic  results  from  sto¬ 
chastic  formulations.  The  computational  limitations  of  Monte  Carlo 
simulation,  especially  high  resolution  or  for  combat  between  large  units, 
are  well  known.  The  work  of  Clark  shows  (see  pp.  102-103  of  reference  12) 
the  complexity  of  an  analytic  solution  for  even  an  idealized  combat  situa¬ 
tion  modelled  as  a  continuous  parameter  Markov  chain  (see  reference  12 
for  further  references).  Moreover,  Bonder  and  Farrell1101  report  excellent 
agreement  between  Monte  Carlo  or  stochastic  simulation  results  and  those 
for  a  corresponding  deterministic  Lanchester-type  model  with  numerical 
solution  generated  by  finite  difference  methods  (see  reference  42).  Of 
course,  verification  of  such  models,  as  with  any  combat  model  is  an  unre¬ 
solved  question  (see  reference  9;  further  references  are  given  in  reference 
38). 

2.  There  is,  moreover,  far  from  universal  agreement  as  to  what  are  the 
significant  variables  that  describe  the  combat  process  and  can  be  used  to 
predict  its  outcome.  For  some  other  views  see  references  19  and  24. 

3.  The  influential  military  philosopher  of  the  19th  century,  Carl  von 
Clausewitz  (1780-1831)  stated  in  his  classic  work  On  War  ( Vom  Kriege), 
“The  best  Strategy  is  always  to  be  very  strong,  first  generally  then  at  the  de¬ 
cisive  point.  .  .  .  There  is  no  more  imperative  and  no  simpler  law  for 
Strategy  than  to  keep  the  forces  concentrated,”  (p.  276  of  reference  13). 

4.  We  note  that  for  any  homogeneous  force  model  of  the  form  dx/dt  — 
—  h(t)A(x,y),  dy/dt=  —h(t)B(x,y),  the  transformation  r=  j0‘  h(s)  ds 
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(where  we  must  place  restrictions  on  h(t)  to  ensure  integrability)  leads  to 
the  autonomous  system  dx/dr=  —A(x,y),  dy/dr=  —B(x,y),  where  A, 
B,  and  h  are  nonnegative  functions  of  their  argument (s). 

5.  Bonder  and  Farrell1101  use  range,  i.e.,  force  separation,  as  the  independ¬ 
ent  variable  in  their  work,  while  Taylor1381  uses  time,  as  we  have  done  here. 
Moreover,  the  X  force-level  equation  (9)  with  attrition-rate  coefficients  (6) 
is  apparently  a  differential  equation  that  has  not  been  encountered  previ¬ 
ously  in  applied  mathematics.  The  differential  equation  under  considera¬ 
tion  could  not  be  found  among  the  445  linear  second-order  equations  tabu¬ 
lated  in  reference  22  or  the  596  tabulated  in  reference  27. 

6.  One  could  equally  well  develop  a  generalization  of  the  constant-coeffi¬ 
cient  solution  expressed  as 

x(t)  =  {[zo+2/o(a/5)’]exp{  -  (ab)H]  +  [x0—yo(a/b)i]exp{  (ab)1^  }/2. 

In  this  case  two  of  the  general  Lanchester  functions  would  be  like  decreasing 
exponential  functions  (with  the  appropriate  arguments)  and  the  other  two 
increasing  exponential  functions.  We  have  chosen  to  generalize  the  form 
of  (3),  since  it  appears  to  be  more  convenient  for  parametric  studies  in 
which  one  might,  for  example,  want  to  vary  initial  force  levels  or  some 
measure  of  relative  effectiveness.  (In  the  constant  coefficient  case  the 
latter  is  a/b.) 

7.  The  singularities  of  the  X  force-level  differential  equation  (9)  occur  at 
the  zeros  and  singularities  of  a(t)  and  at  the  singularities  of  a(t)b(t).  For 
example,  if  a(t)  is  a  rational  function,  i.e.,  the  ratio  of  two  polynomials,  then 
its  logarithmic  derivative,  { 1  /a(t)\da/dt,  has  a  simple  pole  at  the  zeros  and 
poles  of  a(t).  Consequently,  when  a(t)  and  5(f)  are  both  rational  func¬ 
tions,  t0  belongs  to  the  set  of  points  consisting  of  the  zeros  and  the  poles  of 
a(f)  and  5(f).  Finally,  we  note  that  t0=—C  for  (9)  with  the  attrition- 
rate  coefficients  (5)  and  also  (6). 

8.  W (xi,  Xi)  may  vanish,  however,  at  <0,  which  is  a  singularity  of  the  dif¬ 
ferential  equation  (9)  (see  Note7  above  and  pp.  141-142  of  reference  29). 

9.  For  a  discussion  of  algebraic  addition  theorems  see  Chapter  II  of 
reference  18.  Hancock  gives  the  following  theorem  (p.  37,  reference  18): 
a  necessary  and  sufficient  condition  for  a  single-valued  analytic  function 
/(z)  to  possess  an  algebraic  addition  theorem  is  that  there  exist  between  the 
function  f(z)  and  its  first  derivative  an  algebraic  equation  whose  coeffi¬ 
cients  are  independent  of  the  argument  z.  No  more  comprehensive  result 
could  be  found  in  the  more  recent  monograph  by  Aczel.[2i 

10.  A  function  similar  to  F,( x)  was  introduced  by  Ludwig  Schlafli  in 
1867 1321  and  another  appears  in  a  posthumous  fragment  by  William  King- 
don  Clifford  (1845-1879)  (see  pp.  346-348  of  reference  14).  Greenhill1161 
suggested  that  such  a  function  would  be  convenient  for  certain  engineering 
problems,  although  he  was  severely  criticized  by  Watson  (see  p.  91  of  refer- 
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ence  43) .  Since  no  tabulations  of  the  above  or  modified  Bessel  functions  of 
the  first  kind  of  fractional  order  (except  for  the  restrictive  set  of  values 
v=  ±}/i,  ±yi,  ±Ki  =b%>  ±/4)  exist,  it  would  seem  appropriate  to  intro¬ 
duce  the  LCS  functions  as  we  have  done. 
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